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ABSTRACT 



Recently, exponential models of refractivity have been 
devised by several investigators to replace the linear, "standard" 
refractive atmosphere. In this study, average values of N are 
evaluated from the data at the mandatory radiosonde levels up 
to 25 mb^ and mean coefficients satisfying A/ U/ = /V. ) 

have been evaluated for a group of eight mid-latitude stations, 
and for the entire winter season of 1958-59. The coefficients 
plotted as functions of the height of the i-th layer- 
center, possess certain characteristic features in common at all 
stations . 

The writer wishes to express his appreciation for the 
assistance and enoouragement given him by Professor Frank L. 
Martin of the U. S. Naval Postgraduate School in this investiga- 
tion . 
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!• Introduction 



The refractive index n is defined as the ratio of the 
speed of an electromagnetic wave in a vacuum to the wave speed 
in the medium concerned. However, this particular form of the 
refractive index is not vo\y convenient in the iua of 

ray-paths; hence the modified refractive index, N is normally 
used . 

The refractivity N, which is equal to (n-1) x 10 b , may 
be obtained from radiosonde observations using the equation 



N - y ( p * ) = (*-' ) x lo 6 

where T is the temperature in degrees Kelvin, p is the 
pressure in millibars, e is the vapor pressure in mbs, and 
A and & are constants. 

Using the Smith 6c Weintraub [4j values of these 
constants. Equation (1) may be written: 



N 



11. Q P 



313 25 & x 

T * 



( 2 ) 



A Standard* 1 distribution of N in the vertical may be 
obtained by averaging over a sufficient number of individual 
soundings. Bean and Thayer flj have suggested that its 
curve can be fitted by an exponential formula. The model 
which has recently been adopted by the National Bureau of 
Standards is of form: 



~c 2 



N - A/ s e 

where N is the refractive index at any altitude 



(3) 
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N s is the mean surface value of the refracts ’ < r*, ; . 

z is the height above the surface, 

c is an empirical constant, affording ’’best fit” in 
Equation (3). 

Equation (3) is used only with the mean values of N 
averaged over a long enough period to eliminate all transitory 
trapping layers. The study by Bean and Thayer [ij also suggested 
further study by obtaining individual values of N s and c for 
varying climatic regions. This was the objective of the 
first approach of this author to the problem. Later, however, 
it was decided to make a finer-scale determination of the 
parameter c by determining a mean-value cj_ for each layer 
between mandatory levels. Thus a working equation of form 



was employed. Here c^ is temporarily assumed to be constant 
in the i-th layer, which extends from to and N£ is 

the value of N at z = Note finally that c£ may be 
obtained by differentiation of (4) as 




-c: O- 2 \) 






(4) 
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- _1 d N 

A/ d Z 



(5) 
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2. Data Processing 

The stations utilized in this study were selected to 
represent a wide variety of climatic regions. The cj~ values 
were determined from the three-month means of pressure, tempera- 
ture, and moisture content, for all the mandatory levels at 
these stations during the winter of 1958-59. These data were 
obtained from the radiosonde checked data found in the 
Northern Hemisphere Data Tabulations Daily Bulletins [2] , 

The use of data obtained from radiosondes introduced some 
uncertainty of the relative humidity when ’’motorboat ing M was 
reported. However, motorboating usually occurred near 400 mb, 
and there the error introduced by omitting ’’motorboat ing" 
reports from the mois ture-averaging process was found to be 
£ 0.1 mb. The twice-daily observations taken at 0000 and 
1200 Greenwich times were averaged. This served to filter 
out the diurnal variation, and to make the stations, some of 
which were separated by as much as six hours time-difference, 
more compatible. Two methods were used in obtaining the 
values utilized for the determination of the surface refractive 
index data. When the surface report was part of the published 
radiosonde, this was the preferred data-source. When this 
radiosonde level was not available, it was necessary to employ 
the surface observation for the station. However the surface 
observation was recorded only once each day, at 1200Z, and 
this would result in a comparison of a once-daily observation 
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against a twice daily one when going from the surface to 850 
millibars. As a test, the surface three-month mean N-value for 
Washington was computed using both methods and a difference of 
only one N-unit in 300 was noted. Since Washington had one of 
the largest surface diurnal variations of the stati-ons studied, 
this indicated that the once-a-day surface observation reports 
could be substituted for the two-a-day radiosonde surface 
reports without serious error. 

For the actual determination of mean values, the standard 
levels from the surface up to 25 mbs were used. The twice-daily 
observations at these levels were totaled and then divided by 
the total number of entries. All the calculations were done on 
either tape adding machines or hand calculators. A maximum of 
180 observations were used for each lev*], provided ‘that?- nv dtf-ta 
were The mean vapor pressure was obtained by determining 

the vapor pressure from each observation using the temperature 
and relative humidity data^ and finding the mean of these. For 
the majority of the stations, moisture was recorded only up to 
the 400-mb level, at which level the vapor pressure was found to 
be less than 1 mb. The only station that had moisture data for 
the 300-mb level showed the mean vapor pressure to be less than 
one-tenth of a mb. This showed that very little accuracy would 
be lost by not including moisture considerations above 400 mbs. 

After the three-month mean values of pressure, temperature, 
and vapor pressure were obtained for the mandatory levels of 
each station, the three-month mean N-j_ values were determined 
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for these levels using Equation (2). These values and the 
appropriate pressure levels are given in Table 1. 



The -values for each layer were obtained using the 



refractivity equation in the form: 



&sr\ 



N\ 



( 6 ) 



2u,-Zi 

where -z^ is the average-value of the thickness in the 

standard layers employed in this study. 

The values of c^ for each station were then plotted on two 
cycle semi-logarithmic paper versus the midpoint- height of each 
layer as ordinate. These midpoints were taken from the U. S. 
Standard Atmosphere elevations for the stated pressure. A 
discussion of these results and of some applications are given 
in the following sections. 
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Table 1. Three month mean-values of at indicated pressure -levels , and stations employed in this 
study(in N-units)? 
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3* Discussion and Results 



(a) Typical Model 

By Equation (5), C£ is a good indication of the vertical 
distribution of N, An increasing value of ci indicates a 
larger rate of decrease of N with height. Most stations 
studied showed two main values of c^ in the vertical. A 
nearly constant value was indicated for the troposphere, and a 
somewhat larger value shown for the stratosphere. These 
values were separated at most stations by a characteristic 
feature, which will be referred to as a dual discontinuity . 

The dual discontinuity is evidently associated with the tropo- 
pause; it can be described as having a minimum value of c £ 
and immediately above it a pronounced increase to a maximum 
value. A model of cj[ vs height which depicts schematically 
the association between the dual discontinuity and the 
related temperature sounding is shown in Figure 1, which may 
be regarded as typical of most mid-latitude stations. 

In order to understand the association depicted in 
Figure 1, consider N as given by Equation (2) with e^O. Then 



where ft* A is the autoconvec t ive lapse rate, and # is the 




( 7 ) 



Partial differentiation of N in (7) with respect to z, 



and use of the hydrostatic equation, leads to 
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A 




Figure 1, Schematic relationship between the. radiosonde 
and the associated refractivity profile, for mid- 
latitude stations. 



atmospheric lapse rate. By use of Equation (5), it follows 
that 



Hence c^ is relatively large in a thermally stable layer 
and relatively small in a thermally unstable layer. Thus if 
we disregard moisture, any minimum value of Cj[ is associated 
with a relatively steep lapse in the layer, and a maximum 
value with a relatively stable lapse. The discussion of the 
individual stations, which follows, will point out some 
deviations from the typical model, Figure 1. 

One can also qualitatively estimate what the effect of 
moisture lapse would be on c^. Empirically one knows that a 
large negative value of c £ a f/&z may he associated with a 
strong hydrolapse and this would cause c^ to be large 
[^Equation (5)] . On the other hand, if the temperature 

lapse is not abnormal, may be a minimum as a result of a 
small lapse (or possibly an increase with height) of water 
vapor. Moisture effects may be associated with some of the 
deviations of the individual c ^ graphs from the ideal one 
of Figure 1. 

(b) Adak 

The tropospheric distribution of c^ vs height, 

Figure 2, showed a gradual decrease of c ^ with increasing 
height up to 9 km. Above this level, a pronounced dual dis- 
continuity was found, indicating the presence of a tropopause 




( 8 ) 
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layer. A second dual discontinuity of lesser magnitude was 
noted above the first. This was taken as an indication of a 
multiple tropopause, which is in keeping with the model of 
Figure 1. 

(c) Denver 

The graph for Denver, Figure 3, shows a smaller value of 
c i in the troposphere than in the stratosphere. There is a 
sharp increase in c^ to a maximum value near 13 km, with no 
definite dual-discontinuity minimum below. 

(d) Ship "P", Tatoosh Island, Kodiak, Washington, D.C. 

The graphical results for these stations are shown in 

Figures 4, 5, 6, 7, respectively. These graphs have been 
grouped together because they all have features closely re- 
sembling the idealized model of Figure 1. The tropospheric 
C£ decreases slightly with elevation, and is followed by a 
faster rate of decrease ending in the region of 8*10 km. 

The stratospheric c^ then conforms rather closely to the 
schematic model. Kodiak shows a lesser tendency to have the 
pronounced minimum c^ as part of its dual discontinuity. 

This might well be a characteristic of its higher latitude, 
together with a greater variability of tropopause height. 

(e) Ship n N M 

The graphical results for this station, Figure 8, shows 
a greater decrease of tropospheric c^ with elevation than is 
shown in the idealized model. This is due primarily to the 
greater moisture content, and consequently greater N-value, in 
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low levels at this subtropical maritime station. Thus the low- 

level values of are larger than those for most other stations, 

although the minimum just below the tropopause has nearly the 

-5 

same value (about 11.5 10 per meter). 

(f) Lake Charles 

The graph for (pi/o* Charles, Figure 9, shows the tropo- 
spheric ci -value decreasing to a minimum which extends from 

6 km to 10 km. Above this layer^ c^ increases to n wvv/ ly 

* 

constant maximum value in the stratosphere. The - WJ. r.v tant 

c^ in tlft. range 6-10 km is a characteristic of variable tropo- 

• 

pa,; neight, during the observation period (c.f. Kodiak, 

figure 6). 

(g) Composite of all stations 

/ 

The composite graph of all stations, Figure 10, conforms 
closely to the idealized model. The cp -value decreases 
slowly with height in the troposphere, and is separated from 
the larger stratospheric c^ -values by the dual discontinuity. 

A secondary tropopause is also indicated by the less pronounced 
upper dual discontinuity in the neighborhood of 15 km, 
within the stratosphere. 
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4, Other Applications and Conclusions 

By means of Equation (6) a method of solving for c i 
between mandatory radiosonde levels was established. The 
resulting value of c^ was there the mean value for the layer 
of thickness A 2 2 ^ • However c^ turns out to be 

somewhat variable with height. Thus in Figures 2 through 
9, all graphs of c^ vs height have been drawn as continuous 
functions. This means that Equation (4) may be reconsidered 

from the viewpoint that c^ be variable with height in any 

♦ 

layer. In this case, however, the mean value c^ between 
standard levels shall be the same as that computed by 
Equation (6). 

In the field of radar-propagation, one frequently must 
deal with the ray-curvature K = ^ ^ 5 given by 

K = /o' fe +.0^18] c^/3 ( 9 ) 

where /3 is the ray-direction relative to the spherically 
horizontal level cj . Up to several years ago, it was 
customary to use as a standard <^ / '^ z ^a value characteristic 
of the low troposphere, namely d /v r ~ ( oC ) where a 

is the radius of earth. We know now from the work of Bean 
and Thayer and others (including this work) that (^ N /sz)d^ 
is negative but decreases in magnitude. Thus, by Equation 
(5) with c^ and N both decreasing with height in the tropo- 
sphere, be relatively small just below the 

tropopause. 

Consider still another refinement in Equation (4), 
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namely that is empirically given as a function of height 
of the type shown in Figures 2 through 9, inclusive, for the 
localities considered. Then logarithmic differentiation of Equation 



(4) gives c8/v 

N d~Z 




or 



^ ¥ — — /V Cj. f ( .4 34?) 




This expression for dP ^ has a correction laccur, as 
compared with Equation (5), given by the second term on the 
right side of (10). Actually this term can be readily com- 
puted in finite differences* from graphs of c-r vs log, c (Z-2;) 
(for example, see Figures 2 through 9). The correction 
factor may take on some importance in radar -propagation in 
the tropopause layer^ in view of the large increase of 
with elevation occurring there. 

Finally it should be noted that if an equation of the 
Bean and Thayer type, Equation (3), is desired. It may be 
shown that c is a weighted-mean of the individual c^, or 



Equation (11) makes it apparent that the use of a mean 
c for the troposphere is reasonable. However, it suggests 
that one should use a different weighted-mean for the lower 
stratosphere than for the troposphere. In Equation (11), 
Z.^, is the height above the surface, so that the averaging 
must extend down to the earth's surface. Equation (11) may 
be readily solved fox c using arty graph .of ?tbe type shown in 



2^ + 1 c 



• * * + C , A z 



( 11 ) 
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Figures 2 through 9, by selecting the c^-values and weighting 
them by the indicated thicknesses AZ; . However, it should be 
remembered that in most cases only half of the lowest layer , 

A 2 , is shown in these figures. For the stations studied 
in this paper, the N s -values appropriate to the Bean and 
Thayer exponential form, Equation (3), are listed in the first 
row of Table 1. 
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